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ABSTRACT 

We study a simple model of compressible reacting flow. First, we derive 
a dispersion relation for the linearized problem, making a distinction between 
frozen and equilibrium sound speed. Second, we study the stability of the Von 
Neumann-Richtmyer scheme applied to this model. One finds a natural generali- 
zation of the C,F,L. condition. 
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INTRODUCTION 


A simple model — one irreversible reaction — of compressible reactive flow 
is presented in this report. Equations of gas dynamics are augmented by one 
equation for the progress variable of the chemical reaction and the appro- 
priate equations of state. The resulting set of equations is a nonlinear 
hyperbolic system with source terms. However, acoustic waves, i.e., small 
perturbation of a given constant state, will not be described any more by the 
wave equation but by a more general dispersive equation. Only in the so- 
called frozen equilibrium limit will the wave equation be recovered with the 
appropriate sound speed. Such a classical analysis (see [ 8 ] , [ 1 ) for instance) 
is extended to the discrete set of equations obtained by using the Von 
Neumann-Richtmyer scheme (see [9]). Such a scheme is known for its simplicity 
and its robustness to handle strong shocks. The known stability results (see 
[9], [7]) for such a scheme are generalized in our context. Both frozen and 
equilibrium limits appear again, and we make a distinction between implicit 
and explicit schemes used to discretize the reaction equation. 


1 . THE DISPERSION RELATION FOR A REACTIVE FLOW 

We consider 1-D reactive flows involving a single irreversible chemical 
reaction A -*■ B. The flows will be compressible, viscous, but heat conduc- 
tion and chemical diffusion effects are neglected. 

Let us introduce the pressure p, the specific volume v, the internal 
energy e, the velocity u, and the progress variable A of the reaction 
A -*■ B, changing from 0 for no reaction to 1 for complete reaction. Using 
Lagrangian coordinates (a,t), equations of motion are given by: 
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We recall that (a,t) are related to the Eulerian coordinates (x,t), where 
x = x(a,t) gives the position, at time t, of a fluid element that was 
initially at position a. We complete the system (1.1) by prescribing p", 
the internal energy (per unit mass) e, and the reaction rate r; precisely we 
have 


( 1 . 2 ) 


p = p+ u)=p-p 
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V > 0 


where pr is the artificial viscosity 


(1.3) 


e * e(p,v,X) 


(1.4) 


r = r(p,v,X), 


In (1.1), we rewrite the third equation in a nonconservative form as 


(1.5) 
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Using (1.1), (1.3) and (1.5), we also obtain 
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( 1 . 6 ) 
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We first notice (see [3] for instance) that 
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Next, we set 
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The first coefficient will be interpreted (see the Appendix below) as the mass 
speed of sound; the second one is related to the themicity (see [3]). 
Therefore, (1.7) reads 

(1.11) H + 

In order to investigate the stability of the solution of (1.1 ) — ( 1.4), we shall 
"linearize" the equation for X around some local equilibrium state 

"ft 

X * X (p,v) (see the Appendix). We have 
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( 1 . 12 ) 


ax _ x (p,v) - x 
at t 


where x is a "relaxation time," assumed given. Next, we linearize the 

equilibrium state by taking the undisturbed fluid span — denoted by the 
subscript 0 — as reference state 


(1.13) 
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Finally, we obtain for the equation governing X 
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The linearized system associated to (1.1) is obtained classically by writing 
v = vo + v', p = po + P'» etc., in ( 1 . 1 ) , ( 1 . 1 1 ) , and (1.14), where 
Vq,Pq,*»» is a constant reference state. Neglecting second order terms and 
performing the approximation of frozen coefficients, we finally find (dropping 
the ' ) 
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_ _o _ * * 

where y , a , d, t, X , X are now given constants. An equivalent form 

P 0 V 0 

of (1.15), which is more convenient for our analysis, is 


(1.16) 
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The dispersion relation associated with (1.16) is classically obtained (see 
[10] for instance) by representing v, u, p, X in (1.16) as a harmonic 


(1.17) v=Vq exp(wt + ika) ..., etc. 

Substituting (1.17) into (1.16), we obtain a linear system for Vq, Uq, Pq, 1q 
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To have a nontrivial solution, the determinant associated with (1.18) should 
be zero, and this gives the dispersion relation, a relation between a) 
and k. A simple computation gives 


(1.19) 
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We rearrange this equation as 
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Of course the roots u> = a)(k) of (1.20) characterize the stability of 
(1.16) according to the sign of their real parts. In view of the results 
given in the Appendix, the coefficients in (1.20) do have a simple interpreta- 
tion in terms of the frozen and equilibrium mass sound speed a^ and "a , 
We refer to the Appendix for the definition; here we just note that 
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We also introduce a new relaxation times 
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2 

Finally, we rewrite the characteristic equation (1.20) in terms of t, a^, 

and *a^ 
e 
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I Of course, as x 00 (resp. x 0), one recovers the limit of the frozen 

(resp. equilibrium) flow. Both limits are dispersion relations associated to 
a wave operator perturbed by a damping term A , with the appropriate 
coefficient, corresponding to the viscous term in (1,2), 

Due to (1,21), one can easily check, using Routh-Hurwicz criterium (see 
I [5] for instance and the next section) 
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(1.24) 


Re io(k) < 0, 


I 

! 

where to = u>(k) is any root of (1.23), Therefore, the solutions of (1.16) 

I 

decay with time, i.e., the solutions of (1.1) are linearly stable. 


2. VON NEUMANN-RICHTMYER SCHEME FOR REACTIVE FLOWS 

Now we investigate the (linear) stability of the Von Neumann-Richtmyer 
scheme used to discretize (1.1) or (1.16) and completed by a simple Euler 
scheme for the X-equation. Let At,Aa denote the mesh scale in time and 
space. For (1.1), we have, using standard convention 




where r n+ ® = r(e P n+ } + (1-0 )p n ...), 0 < 0 < 1. Note that we used 

j + 7 J + 7 j+y 

the energy equation under the form (1.5). Second order accuracy corresponds 
to the choice 0 = ^ , 

Such a scheme is now applied to the linearized version of (1.1), that is 
(1.16), again using the assumption of frozen coefficients. To avoid 

fractional indexes, we make the shift i+^. i, n +^. -► n+1, etc. 

We find 
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where 


p^+® = 0p? + * + (1-8 )p?, etc., 0 < 9 < 1. To carry out a stability 
J 3 3 


analysis for the scheme (2.2), we set 


(2.3) 
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and similar relation for u^, p?, X?. Substituting (2.3) into (2.2), we 
find a set of linear equations for Vg, Ug, Pg, lg* 
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The characteristic equation associated to this system is obtained by computing 
its determinant. We find, for g = sin 2 ^2- 
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where we used the definition of the frozen and equilibrium mass sound speeds 
(see the Appendix) and the definition (1,22) of t\ 

Equation (2.5) is to compare to its continuous analog (1.23). Of course, 
r-1 is one trivial root. To insure stability, we have to check that the 
remaining roots of (2.5) do have modulus smaller or equal than 1. This will 
be achieved assuming conditions on At and Aa. 

To simplify notations, let us introduce 
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Therefore (2.5), up to a factor .Eli , can be rewritten as 


(r-1) [(r-1) + A(r-l) + Br] 


(2.7) 
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We are now looking for conditions insuring that the roots r of (2.7) 
satisfy |r| <_ 1. Classically, we perform the transformation r = 

( | r | jC 1 C ■ > Re z < 0) and apply the Routh-Hurwicz criterium. From 

(2.7) , we deduce 

3 

(2.8) l a.z-’ = 0. 

j=0 J 

The coefficients in (2.8) are given by 

(2.9) a 3 = 2(4 - 2A - B) + (28 - 1)C(4 - 2A - B ' ) 

(2.10) o 2 = 4A + 2(28 - 1)AC + C(4 - 2A - B') 

(2.11) Uj = 2B + (20 - 1)B'C + 2AC 

(2.12) a Q = B'C. 

The Routh-Hurwicz criterium insures Re z < 0, for any root z of (2.8), as 
soon as (see [5]) 


(2.13) a 2 , 0 

(2.14) a l a 2~ a 0 a 3 — 

By the definitions (2.6) (we recall |5 = sin^ and the property (1.21), 


we have 
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(2.15) A, B, B', C > 0 and B > B'. 

To study (2.13) and (2.14), we discuss according to the sign of 20-1 
(0 _< 9 <_ 1). 

i) The case 20 - 1 _> 0. 

Here the scheme for the A-equation is unconditionally stable. In view 
of (2.9) ,(2.10) , and (2.15), the condition (2.13) is satisfied as soon as 

(2.16) 4 - 2A - B > 0 and 4 - 2A - B' > 0. 

Due to (2.6) and (1.21), we finally obtain the sufficient condition 

(2.17) (¥ f t4) 2 + 2 » < I* 

Thus, we recover the usual condition of stability for the Von Neumann- 
Richtmyer scheme (see [9]) but with the appropriate sound speed. Next, we 
look at the condition (2.14), which reads 

(4A + 2(20 - 1) AC)(2B + (20 - 1)B'C + 2AC) 

(2.18) 

+ 2AC 2 (4 - 2A - B') + 2(4 - 2A)(B - B')C > 0. 

Since 20 - 1 0, (2.18) is satisfied as soon as (2.17); hence, (2.16) 

holds. 

In summary, (2.17) insures (linear) stability for (2.2) if 9 <_ 1 . 
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ii) The case 20 - 1 < 0. 

Here the scheme for the X-equation becomes explicit. Again using (2.6), 
we rewrite (2.13) as 

(2.19) 2(4 - 2A - B - 4(1 - 20)C) + (1 - 20)C(2A + B') > 0 

(2.20) 2A(2 - (1 - 20)C) + C(4 - 2A - B") > 0. 

Thus, due to (2.15), (2.19) is satisfied as soon as 

(2.21) 4 - 2A - B - 4(1 - 20)C > 0. 

With (2.6), we thus obtain the sufficient condition 

(2.22) (a f s| )2 + + (1 “ 20)(^-) 2 ~< 1; 

(Aa) a e x 

that is the condition (2*17) enforced by an extra term coming from the 
explicit scheme used for the X-equation, We remark that this extra term 
vanishes in the frozen limit T + + 

Of course, (2.21) implies (2.16), since B^B^; therefore, (2.20) is 
satisfied as soon as 

(2.23) 2 - (1 - 20 )C > 0, 


and we find the natural condition 
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(2.24) (V-L^Ut < 

a 2 t 
e 

It remains to check (2.14), i.e., (2.18). Due to (2.23), it is sufficient to 
satisfy 

(2.25) 2B - (1 - 20 )B'C > 0 
or 

(2.26) At < 1. 

2t 

o o 

But, since a f > a^, (2.26) is a consequence of (2.24). 

In summary, (2.22) and (2.24) insure (linear) stability for (2.2), if 

0 £ 9 i 7 • 

The above results could be easily generalized to similar 2-D reactive 
flows discretized by the natural extension of the Von Neumann-Richtmyer scheme 
(see [6], for instance). Similar results are anticipated for the Godunov 
scheme, which reduces, after linearization, to the Courant-Isaacson-Rees 
scheme ([7]). They will be published elsewhere. 
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APPENDIX 

Here, for the convenience of the reader, we recall some classical results 
concerning the definition of the sound speed in reactive flows (see for 
instance [ 3] , [8] , [ 1] ) . Here we take "p = p, that is, (see (1.2)); we 

suppose y = y = 0. 

According to the laws of thermodynamics. 


(A.l) TdS = de + pdv - vdX, 

where T, S, v are the temperature, the entropy, and the chemical poten- 
tial. From (A.l) and (1.3), we deduce 

(A. 2) Tds = dp + (|^ + p)dv + “ v)dX. 

Now, we consider two cases. 


i) The flow is frozen : dS = dX = 0. 

Then, (A. 2) becomes 
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Therefore, the definition (1.9) of a 1 2 (where ? = p) is that of 

2 2 

frozen mass sound speed p c^. We shall use the standard notation 

2 _ 2 2 

a f = p v 


the 


ii) The flow is in chemical equilibrium : dS = 0 and A = A*(p,v). 


In particular v 0 and dA - 7^-- dp + ■— dv in (A. 2); collecting 
terms we find 
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The equilibrium mass sound speed is then a 2 = pc 2 , and we easily check 
(see (1.20), (1.8), (1.10)) that 
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Similarly, we find (see (1.20)) that 
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(A.8) 
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Finally, in view of (A. 4) and (A. 6), we have 
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